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ABSTRACT: We have developed a matrix calibration procedure that uniquely relates the electric 
fields measured at the aircraft with the external vector electric field and net aircraft charge. Our 
calibration method is being used with all of our aircraft/electric field sensing combinations and 
can be generalized to any reasonable combination of electric field measurements and aircraft. We 
determine a calibration matrix that represents the individual instrument responses to the external 
electric field. The aircraft geometry and configuration of field mills (FMs) uniquely define the 
matrix. The matrix can then be inverted to determine the external electric field and net aircraft 
charge from the FM outputs. A distinct advantage of the method is that if one or more FMs need 
to be eliminated or de-emphasized (for example, due to a malfunction), it is a simple matter to 
reinvert the matrix without the malfunctioning FMs. To demonstrate our calibration technique, 
we present data from several of our aircraft programs (ER-2, DC-8, Altus, Citation). 
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1. Introduction 


Measurements of electric fields within clouds have been made by aircraft for many years 
[e.g., Gunn et at, 1946; Blakeslee et al., 1989; Winn, 1993]. Retrieving electric field 
components from the raw aircraft FM data is perhaps the most difficult aspect of the 
measurement process [e.g., Jones, 1990; Koshak et al., 1994; Koshak, 2006a; Koshak et al., 
2006b]. The electric field as measured by an instrument on the aircraft includes linear 
components from the external electric field (e x ,e Y ,e z ), charge on the aircraft ( e Q ), mill DC 

offset, and various other processes, including measurement errors and non-linear terms (/). 
Typically we use electric field mills (i.e., mills) [Bateman et al., 2006] as our electric field 
measuring instruments on aircraft. The equation relating the output of the i mill (a/) to the 
external electric field components can be represented by: 

a i = k^ix e x ■*" M iY e Y + M jz e z + M^q + M iD + y , (1) 

where the various M coefficients are the responses of the mills to the external field and charge 
on the aircraft, M io is the DC offset of the z' th mill, and we have neglected mill measurement 
errors. 

For most techniques [e.g., Koshak et al, 1994; Koshak, 2006a; Koshak et al., 2006b], 
including ours, we neglect the non-linear terms so that the system can be solved with a linear 
algebraic approach. We also will neglect measurement errors in this analysis. That is, neglecting 
y results in a linear set of equations between the mill outputs and the external electric fields. If 
you expand (1) in time and for all mills on an aircraft, the resultant system of equations are: 
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( 2 ) 


#1 (/) — M xx e x (t) + M lY e Y (0 + z e z (0 -^i q^q (0 ^id 

a 2 (0 = M 2 x e x (0 ■*" M 2Y e Y (/) + M 2Z e z (t) + M 2 q€q (t) + M 2D 

a m (0 = M mX e x (0 + M mY e Y (0 + ^ mZ e Z (0 + ^ mQ e Q (0 + ^ mD > 

where wz is the number of mills on the aircraft. 

The mill DC offset are due to a number of factors that can change from day to day and 
even during a flight, but are usually small and are not significant during storm overflights. 
However, the offsets can be significant when compared to fair weather fields (when we usually 
do our calibrations). The reason why we explicitly call out this term will become clear later in 
the description of our calibration procedure. 

The set of equations in (2) can be represented by a matrix equation: 

a(0 = Me(f), (3) 

where a(i) is the vector of mill outputs as a function of time: 

a, it) 

a 2 (t) 

a(0 = a 3 (0 , (4) 

a m (t) 

e(i) is the vector external electric field (including the charge on the aircraft) as a function of time: 

e X^f) 

e Y (t) 

= (0 ’ (3) 

e Q (t) 

1 

and M is the m x 5 calibration matrix: 
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( 6 ) 



M lx 

M„ 

M lz 

m iq 

M w 

M = 

M 2X 

m 2Y 

M 1Z 

M 2Q 

M 2D 


MnOC 

M m y 

M mZ 

M m Q 

M mD 


Note the last row in e(7) and the last column in M. As discussed above, these factors are 
introduced to account for (and thereby eliminate the effect of) mill offsets in our calibration 
procedure. How these terms affect the calibration will be discussed in later sections. For clarity, 
note that we use capital letters for matrices and small letters for vectors throughout our 
description of our method. 

As given in (3), the matrix M relates the mill outputs to the external electric field and 
charge on the aircraft. The matrix M is constant and unique for each aircraft and mill 
distribution. Given M , one can directly determine the outputs of the mills on the aircraft based 
on the external electric field vector. However, what we have are the mill outputs, a (t) , and what 
we want is the external electric field vector, e(t ) . Pre-multiplying (3) by the transpose matrix 

M r and pre-multiplying again by the square matrix (M 7 M) gives: 

[(M r M) _1 M r ] Me(0 = [(M r M) _I M r ] a (t) 

Ie(0 = (M r M) _1 M r a(0 (7) 

e(0 = Ca(0 , 

where I is the identity matrix and Cs(M r M)' 1 M r which is the Moore Penrose (MP) 
‘pseudoinverse’ [Penrose, 1955] of M. In Matlab®, the function call is C = pinv( M). All of 
our mathematical analysis was done using the Matlab® software package. For simplicity, we 
will use the notation one would use in the Matlab® command line whenever possible. However, 
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we will expand the very compact notation when needed to explain what mathematical 
manipulations are actually used by the software package. 

The elements of C are similar to the elements of M in that they relate the mill outputs to 
the external electric field (with net aircraft charge). For example, e x (t) is related to the mill 
outputs by: 

e X (0 = (0 Cx 2 a 2 (0 + C X 3 a 3 (t) ■+ t- C Xm a m ( t ) , (8) 

where the various C’s are the first row of the C matrix. The electric field solution in (7) is 
termed the least-squares solution since it is the minimum of the scalar s (e) = (a - Me) 2 . 

2. Calibration Procedure 

We use several techniques to minimize the errors in the retrieved elements of M , which 
we will describe here. We have used this procedure to calibrate the FM systems installed on 
various aircraft (e.g., ER-2, DC-8, Citation, and Altus). The general steps are illustrated in Figure 
1. Our overall strategy is to compute an initial estimate of M given by M 0 , and then perform 

j = l,...,P subsequent iterations to obtain refined estimates of M (i. e., M,,M 2 ,...,M P ), where 
M p is our best estimate of M . 

To compute M 0 , we first mount the mills on the aircraft in appropriate locations (step 1 in 
Figure 1). The next step is to fly the aircraft through a “known” electric field (such as the fair 
weather atmospheric field) while putting the aircraft through a series of roll and pitch calibration 
maneuvers and recording mill outputs (step 2). The next step is to make an ideal estimate of the 
external electric field in the aircraft frame-of-reference (step 3 in Figure 1) based on the model 
field profile (fair weather or other know field conditions) and the aircraft roll and pitch 
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maneuvers. The ideal electric field in the aircraft frame is combined with the mill outputs to 
compute M o (step 4 in Figure 1). 

Next, M o is corrected for any known matrix symmetries (step 5). The corrected M o 

matrix is inverted and combined with the mill outputs to create an estimate of the true electric 
field (step 6). This est im ate of the electric field is corrected for any known problems (step 7 in 
Figure 1). It is then used to create the next iterate M, matrix (step 4 of Figure 1). The j = j + 1 

iteration loop in Figure 1 is repeated until M . converges. The M y . matrix usually stabilizes after 
about five iterations, but the general convergence conditions have not been proven or guaranteed. 
We then make an absolute calibration of the matrix to create our best estimate of the true M 

matrix. The details of each of these steps will be covered in the following sections. 

2.1. Mill Locations (Step 1) 

Obviously, the very first step in the calibration procedure is to choose good locations for 
the field mills (step 1 in Figure 1). A proper placement of mills on an aircraft leads to an accurate 
calibration, while poor placement leads to excessive calibration errors which can totally mask the 
true M matrix. In our method, each mill should measure at least one component of electric field, 
and each component of the field should be dominate in at least one mill. Ideally, the mills should 
be distributed so that no two mills have similar sets and polarities of dominate field components 
and, if possible, both polarities of each component should be represented. That way, each mill 
will contribute to the measurement of the external fields. 

An example mill placement for an ideal aircraft is shown in Figure 2 (top example ‘A’). 
Note that each mill has at least one dominate component of the external electric field (and net 
aircraft charge) and each component is dominate in at least one mill. Further note that different 
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polarities are present for many of the field components. The polarity changes make it easier to 
remove the effects of aircraft charge. This placement of mills would make the calibration process 
easier. An example of a poor mill placement is also shown in Figure 2 (bottom example ‘B’). 
Although the mills have dominate components in e x , e Y and e z , they all have approximately the 
same dominance in e Y and e z . Since all mills would then react similarly to e Y and e z , it would 
be very difficult to determine the e Y and e z fields from this mill placement. 

2.2. Calibration Data (Step 2) 

Although any set of measurements of external electric fields can be used for calibration, it 
is best to pick data when the field values can be easily estimated. For example in our 
calibrations, we usually use field data during periods when the aircraft is performing roll and 
pitch maneuvers in fair weather conditions. Under such conditions, the external vector electric 
field is predominantly vertical with a profile that can be estimated \e.g., Gish, 1944]. The roll 
and pitch maneuvers map the assumed vertical fair weather field into the aircraft frame e x ,e Y 

and e z components. An example of a set of calibration maneuvers for one of our aircraft is 
shown in Figure 3. Figure 4 shows the mill outputs during the roll and pitch maneuvers in Figure 
3. 

Our calibration technique works best when the various field components in the aircraft 
frame of reference are distinctly different and varying. Such times are when the aircraft is 
actually performing the pitch and roll maneuvers. Therefore, to maximize the convergence of our 
technique, we attempt to use field values only when they are changing due to the pitch and roll 
maneuvers. Attempting to fit the ideal fields to the mill outputs during times when the ideal 
fields are constant will often result in the solution fitting the field changes that are not traceable to 
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the roll and pitch maneuvers (such as fields due to charge pockets or thermals). By focusing on 
the roll and pitch maneuvers, we force the least-squares solution to model those waveforms and 
not waveforms from other sources. For example, Figure 5 shows the ideal fields during the 
calibration maneuvers from Figures 3 and 4. The variations in the fields should be readily 
identifiable. A poor choice of calibration data is shown in Figure 6. During this time period, the 
Citation aircraft was not making any significant calibration maneuvers. Note the various features 
in the mill outputs that are not reflected in the ideal fields. Using this time period to calibrate our 
aircraft mill configuration would produce a poor estimate of M very different from the one based 
on the data in Figure 3. 

2.3. First Estimate of the External Electric Field (Step 3) 

Given fair weather fields (derived from Gish [1944]), and the aircraft roll, pitch, and 
altitude information, we can derive first guess or ideal external electric fields in the aircraft frame 
of reference. These are shown in Figure 5 (plots g x , g Y , and g z ) for the calibration example 

used in Figures 3 and 4. To determine the g Q coefficient, we need to have charge induced on the 

aircraft frame. Often this is done with a high voltage stinger probe [e.g., Vonnegut et al., 1961], 
however, charges induced on the aircraft due to the maneuvers or power changes can be used if a 
stinger is not available. In addition, we have found that the fields at the mills are often 
dominated by the charge on the aircraft even while flying in fair weather conditions. Since the 
other components are small and are often alternating polarities in different mills (some mills have 
positive components, some have negative), the sum of the mill outputs is often almost exclusively 
due to the charge on the aircraft. Therefore, we can use the average mill value (summed over all 
mills) as a first guess to the fields at the mills due to charge on the aircraft ( g Q ). Note that there 
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is no “absolute” value of g Q . All that is needed is a relative value so that the effects of charge on 
the aircraft can be separated from the other components. The first guess as to the field due to 
charge on the aircraft are also shown in Figure 5 (plot g Q ). 

2.4. The First Estimate of M (Step 4 with j = 0) 

To produce the first guess of the M matrix (M 0 ), we gather the mill outputs as a function 
of time into an mxn matrix A : 

A = [a(f,) a (t 2 ) a (t 3 ) ••• a(t„)] , (9) 

where each a(/-) is the column vector of m mill measurements at time t- and n is the number of 
observation times used from the calibration period. We then collect the first guess of the external 
fields (including aircraft charge) into a 5 x n matrix G : 

G = [g(0 g (* 2 ) g(f 3 ) g(U] • (10) 

Each column vector g(/) is a first guess of the electric field in (5) at a particular time, and is based 
on the field profile and aircraft roll and pitch maneuvers, and n is as described for (9). It can also 
be called the "ideal" fields as G is the field that would be present if the field profile was correct 
and there were no other contributions to the electric field (such as charge pockets). We then 
obtain M 0 by Matlab® right matrix division: 

M c = A/G , (11) 

which is a short-hand representation of the least squares solution [Anderson et al., 1999] for 
finding M 0 . The definition of this matrix division is provided in the Appendix. 


10 



We have found that the most common source of non-field components are DC offsets in 
the mill outputs. To prevent the mill DC offsets from distorting our calibration, note that the fifth 
row of G contains all Is, and the contributions to the mill outputs from the offsets (which are 
constant over the time of the calibration) are hence contained in the fifth column of M 0 . DC 
offsets cannot easily distort the minimization of the other components because of this 
arrangement of constant terms. This approach is used for all subsequent iterations to prevent the 
mill offsets from biasing solution M p . 

2.5. The / h Estimate of E (Steps 5 and 6) 

Once we have the M 0 matrix (including the offset terms), or more generally any iterate 
M . , we correct it for any know aircraft symmetries (step 5). For example, a mill mounted along 
the centerline of an aircraft should not have significant e Y components. We call the corrected 
matrix M* . In step 6, we take the MP pseudoinverse to obtain: 

Cj = pinv(M*j) ; j = 0,l,...,P. (12) 

Invoking (7), we estimate the true external electric fields with: 

F,.=C,A , (13) 

where F . is the 5 x n matrix of the j th estimate of the true external electric fields: 


A 


Aft) 

Aft) 

Aft) 

Aft) 

Aft) 


Aft) Aft) ••• Aft) 
Aft) Aft) ••• Aft) 
Aft) Aft) Aft) 
Aft) Aft) •" Aft) 
Aft) Aft) ••• Aft) 


(14) 
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The first four rows of F y are the various field components (including aircraft charge) while the 
last row approximates a row of Is (i.e., as M ; . approaches M, the last row of F, approaches a 
row of Is). 

Although the ideal electric field, G , is a good first guess, the estimated electric field, F. , 
often have temporal and spatial variations that are not found in the ideal fields. For example, 
Figure 7 compares G (ideal fields, gray waveform) with F 0 (estimated fields, black waveform) 
during a calibration run. Note that the fields are much more variable in the “true” case. Many of 
the “bumps” in F 0 correspond to roll and pitch maneuvers (as indicated in the G waveform), 

however, there are a number of excursions in F 0 that are not caused by the roll and pitch 
maneuvers. Attempting to force the M y matrix to fit G when F ) is more like the true fields can 
introduce significant errors in the matrix determination process. 

2.6. Correction of F y (Step 7) 

The next step is to correct the components of F y . To do so, we must introduce the true 
electric fields at times t x to t n (the actual e x , e Y , e z , and e Q fields at the measurement times): 


e x (A ) 

e x (A ) 

e x (A ) 

e x (A ) 

e Y (tj ) 

e Y (A ) 

e Y ( t 3 ) 

e Y (A ) 

e z (A ) 

e Z (A ) 

e z (A ) 

e z (A ) 

e Q (A ) 

e Q (A 2 ) 

S Q (A ) 

• e e(A) 

1 

1 

1 

1 


(15) 


which is a 5 xn matrix. Combining (3), (9), and (13) gives: 
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( 16 ) 


F.=C,A = C y [a(0 a (t 2 ) a (f 3 ) - *(/„)] 

= C y [Me(0 Me(/ 2 ) Me(g ••• Me(g] 
= CME 
.. F y = H.E, 


where the 5 x 5 matrix H y is a constant defined by: 


Hz* 



H xq 

H w 

Hy* 

Hyy 

Hy, 

H y q 

H m 

tizx 

H z y 

H zz 

h zq 

h zd 

H qx 

Hqy 

H qz 

h qq 

h qd 

H dx 


h dz 

H dq 

h dd 


(17) 


From the result in (16), one can write any estimated field component, for example, the X - 
component, as: 

f x (ti) = H xx e x (t i ) + H X ye Y (t i ') + H xz e z (t i ) + H X geg(t i ') + , (18) 


where it is understood that this is associated with the y'th iteration; i.e., we drop the j index on the 
/ and H scalars for brevity. 

Note that as M y approaches M, C y in (12) approaches (M 7 M)' ‘M 7 , and hence 
H y = CyM = (M r M) _1 M r M = I . In other words, if M y approaches M, f x {t t ) approaches 
e x {t i ) because the diagonal elements of H y approach unity while the off diagonal elements of 
Hy approach zero. The meaning of (18) is that each component of F y is dominated by the same 
component of E (i.e., the diagonal elements of Hy ), with smaller contributions from the other 
components of E (i.e., the non-diagonal elements of H ; ). By looking for features in the various 
Fy waveforms that are due to the various E components we expect from the calibration 
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maneuvers, we can estimate the contamination of the various components of F ; . , (the non- 
diagonal elements of H . ) and subtract them to produce an F ; . closer to E. 

An example of a F 0 (based on the first guess M 0 and the actual mill outputs) along with 

G fields is shown in Figure 8. We know that the details of E will be different than G, but the 
major transitions will be due to the roll and pitch maneuvers, which will be very similar in both E 
and G. In this example, f x (with j = 0, i.e., the fields based on the first guess M 0 ) has e Y and/or 

e z contamination while f Y and f Q each have e x contamination. We can identify these cross- 
contaminations because we are monitoring the aircraft electric fields during times when we 
approximately know what the fields should be doing (i.e., the G fields). For example in Figure 8, 
during the period 16:14 to 16:15, the aircraft is making pitch maneuvers. We know that during 
that time, the e x and e z fields should be varying while the e Y and e Q fields should be more or 

less constant. Therefore, the “bumps” in the f Q and f Y fields during that time that resemble the 
e x excursions are likely due to contaminations of those components by e x . By comparing the 
amplitude of the “bumps” in the f Y (due to e x contamination), we can estimate the amount of 
contamination and subtract a portion of f x from f Y to produce the improvement f Y (as show in 
Figure 9): 

f; = fr-ccf x , (19) 

where a is an empirically determined factor that is set to best eliminate the cross-contamination. 
As shown in Figure 9, we perform the same operations on the f x , f z , and f Q components to 

obtain the corrected F . , which we write as Fj . 
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One of the reasons this step is done by hand is it is difficult to identify the “bumps” 
associated with the calibration maneuvers under “real world” conditions (as seen in Figures 8 and 
9). Many of the “bumps” (such as the ones in f z and f Q around 16:17) are not readily 

identifiable as due to calibration maneuvers. We have found that these bumps not due to 
calibration rolls or pitches make automatic processing of the calibration data impossible. 

2.7. Subsequent Estimates of M (step 4 with j > 1) 

Next, note in Figure 1 that j is incremented. This means that the F' discussed above is 
now written as F'._, . Hence, the next iterate (and all subsequent iterates) are written as: 

m. = a/f;. i; j=\,.,.,p , (20) 

where A is again the mill outputs given in (9). The process is repeated until the j iterate M , 
converges and the field outputs (F ; .) stabilize to reasonable waveforms (like those in Figure 9). 
The (relative) calibration matrix obtained is M p . Convergence here simply means that the last 
few iterates (M^jM^^M,,) differ by a sufficiently small amount. For typical combinations of 
FM, data, and aircraft, this process usually takes about five iterations, with the majority of the 
iterations used to fully identify the various components in the true electric field (i.e., which 
‘bumps’ in the waveforms are caused by which components of the true electric field). 

2.8. Absolute Calibration (Step 8) 

The steps in the previous sections only provide a relative calibration, that is, the 
components in F p are correctly proportioned to the components in E. However, each component 
may differ (in a unique way) by a multiplicative constant from their true magnitude. Essentially, 
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we have created a F with the correct waveshapes, but not the correct amplitudes. We can 
determine the relative multiplicative constants between the various components (e.g .,f x , f Y , & 

f z ) by checking the relative responses of the various components to the roll and pitch 

maneuvers. However, even after this step, we need one more constant to create a final calibration 
matrix. 

To determine the final multiplicative constant, we need to compare our fields calculated 
with the final M p matrix to a known electric field. This is typically accomplished by flying the 
aircraft near a calibrated ground based mill and comparing the relative amplitudes of the fields at 
the ground and at the aircraft. The ratio of the two fields produces the factors needed in the final 
absolute calibration of the M matrix. Figure 10 shows an example ground-based mill overflight 
for one of our aircraft platforms. 

3. Results 

We have used this calibration method to create M p matrices for four different aircrafts. 

Figure 1 1 show the electric fields based on the final calibration of the FMs on the NASA ER-2 
[Hood et al., 2006]. Figure 12 shows the electric field vector based on the final calibration of the 
mills on the North Dakota Citation [Dye et al., 2004], Figure 13 shows the vector electric field 
components based on the final M p matrix for an overflight of a small storm with the NASA DC- 
8. Figure 14 shows the fields during an overpass of a thunderstorm with the Altus aircraft [Mach 
et al., 2005] based on the final M p matrix for that aircraft and FM combination. In each case, 
the aircraft flew over and slightly to the side of the electrified cloud. 

For such overflights of typical storms (positive charge near the top of the storms), the e x 
field component should first go negative as the aircraft approaches, pass through zero at the 
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closest approach to the storm, reach a positive peak, and then decay away to zero as the aircraft 
recedes from the storm. If the aircraft passes directly over the center of the charge, the e Y field 

should be close to zero. However, if the aircraft passes to the side of the storm, the field trace 
should have somewhat of a Gaussian shape with the polarity depending on which side the aircraft 
passed. The e z field should also have a Gaussian shape with positive polarity. The shape of the 

e Q component depends on so many different factors that it is difficult to characterize. Any 

lightning will be superimposed on these general shapes. 

The wave shapes in Figures 11, 12, 13, and 14 follow these general trends, indicating that 
the calibrations based on our M matrix method are reasonable. Unfortunately, further 
verification of the calibrations requires knowledge of the electric fields that we simply do not 
have. 

4. Mill Output Weighting 

Once the final matrix is calculated, the inversion process can be used to create a C 

matrix with user defined properties. For example, if a C matrix is formed using the following 
formula: 

C = pmv(WM p )W , (21) 

where W is a m x m diagonal weighting matrix. The diagonal elements of W (i.e., 
W H , i = 1, ..., m ) will determine how much the i‘ h mill is emphasized in the resultant C matrix. 
For example, given an aircraft instrumented with 6 mills, we can take mill #3 totally out of the C 
matrix by simply nullifying element W 33 . With the remaining five mills equally weighted, W is 
written as: 
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1 0 0 0 0 0 
0 1 0 0 0 0 
0 0 0 0 0 0 
0 0 0 1 0 0 ’ 
0 0 0 0 1 0 
0 0 0 0 0 1 


( 22 ) 


where the element corresponding to the third row, third column is zero. As long as the original 
M p was correct, the resultant C matrix will create valid fields, but have no contributions from 

mill #3. The diagonal elements of W need not take on binary values (shutting on or off mills), 
but can assume arbitrary values to emphasize or de-emphasize mills in any way pleased based on 
the specific quality of individual mill outputs. For example, to lower the influence of mill #3 by a 
factor of ten, the zero in W 33 in (22) would be replaced by 0. 1 . 


4.1. Example Use of W 

Figure 15 shows the mill outputs during a time when one of the mills was malfunctioning 
(mill QL). Figure 16 shows the field values based on the C matrix with all mills (including the 
malfunctioning mill). Note the significant errors in the plots. All of the field components are 
essentially masked by the mill errors. If we produce a W matrix with element W 22 equal to 0.1 
and use (21) to compute C , we get the data in Figure 17. Note that the influence of mill #2 is 
much less. If we produce a W matrix with element W 22 equal to 0, we get the data in Figure 18. 
Note that the noise from mill #2 is totally eliminated. 

Each C matrix was based on the same matrix but with different W matrices. 
Figure 19 shows the expanded f Q component for each of the three different cases (no weighting, 
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mill #2 reduced by a factor of 10, mill #2 totally removed). Note that the overall field component 
shape is the same in each case. The only difference is the inf luence of the noise from mill #2. 

Reasons for reducing the influence of a mill, instead of simply eliminating it, include 
cases where the mill is critical for determining one component of the external electric field. For 
example, one of our aircraft had four mills mounted in a similar way to mills two through six in 
Figure 2. The aft mill was critical in the estimation of e x , but due to the mounting location, its 
output was often very noisy. By setting its weighting function to 0.1, we were able to 
significantly reduce its influence on the e Y , e z , and e Q components (where it was not needed) 

while still using it to deter m i n e e x (where it was critical). 

4.2. Other Uses of M Matrix Method 

If the real-time analysis program stores M p and calculates C on-the-fly, mill health 
information can be used to modify which mills are used in the determination of the field values 
by manipulating the W matrix in (21). This would allow for real-time changes in the calculation 
of field values based on mill metadata. Individual mills could be turned off and on, or their 
influence simply reduced based on their health in real-time. 

5. Conclusion 

We have developed a method of determining the relationship between the external electric 
fields and the outputs of electric field meters for a generalized aircraft platform. As long as the 
FMs are not placed in a pathological fashion and there are more FMs than electric field 
components, the method converges after a few iterations. By determining the M matrix (the 
matrix that produces the FM outputs from the field values) instead of the C matrix (the matrix 
that produces the field values from the FM outputs), we can use the properties of the M matrix to 
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simplify the calibration process and also add the ability to emphasize or eliminate FM outputs on 
the fly. 
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Appendix 


Solution Forms for the M Matrix 


In this appendix, two standard solution forms for the M matrix are provided. One form 
was given in Koshak (2006a), and the other form is used in this paper an detailed here. It is then 
shown that the two forms are equivalent. This helps better establish the inter-relationship between 
Koshak (2006a) and the current writing. 

In Koshak (2006a) the following solution form for M is provided. For j = aircraft 
orientations of an aircraft having calibration matrix M, one can write the mill outputs a ; = a(? •) 
at time t } due to ambient electric field e y = e(t - ) as 

2Lj = Me . . (A.1) 

Next, each side of this equation is post-multiplied by the transpose , and then summed for all 
j = orientations. This gives: 




7=1 


7=1 


7=1 


M=£ 


T 

a e 

7 7 


7=1 


s 

W=‘ 


T 

e e . 

7 7 


(A.23) 


This last expression in (A.2) is the first expression in (4) of Koshak (2006a), but here we omit 
mill errors and factors of 1 / n for brevity. 
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In the main text of this paper, we have used a different notation. First, we defined the 
following two matrices: 

A = [a 1 a 2 ... a n ], E = [e, e 2 ...e„] (A3) 

which allow us to write all j = 1 ,...,« of the equations in (A.l) as 

A = ME . (A.4) 

Postmultiplying by the transpose E r and inverting gives 

M = AE r (EE 7 ’) -1 . (A.5) 

However, it can be shown by straightforward matrix multiplication that 

AE r =i>. e J , EE T =Y,e/j- (A.6) 

j = 1 7=1 

Therefore, the solution for M in (A.2) given in Koshak (2006a) is identical to the solution in 
(A.5). In addition, note that the right matrix division defined by Matlab® and used in the main 
text (equation 1 1) is defined by (A.5), that is 

M = A/E = AE r (EE 7 ’) -1 . (A.7) 
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Figure 1. Steps used to estimate the M matrix. The process of refining the M ; matrix usually 
takes about five iterations (j — 0, 1, ..., 4) to settle to the point where further changes in the M ; 
matrix produce no discemable changes in the final electric field values. We do not correct the 
My matrix for any know symmetries on the aircraft the final time because the aircraft is not 


exactly symmetric. The best estimate is called M p in the text. The final absolute calibration 
creates our best estimate of M. 
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Side View 


Figure 2. Example placements of mills on an ideal aircraft. Version “A” shows a “good” 
placement of mills. Mill #1 output is dominated by e x and - e z . Mill #2 output is dominated by 
e Y and e z . Mill #3 output is dominated by e Y and - e z . Mill #4 output is dominated by - e Y and 
- e z . Mill #5 output is dominated by - e Y and e z . Mill #6 is dominated by - e x and e z . All mills 
have a strong e Q component. Note that the mills have different dominate components of E. This 
configuration should converge using our calibration method. Version “B” shows a “poor” 
placement of mills. In this version all mills have similar values of e Y and e z , making 
convergence of our calibration method difficult if not impossible. 
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16:13 16:14 16:15 16:16 16:17 16:18 16:19 

Figure 3. Example roll/pitch maneuver set obtained during the Altus aircraft calibration flight. 
The roll and pitch angles are in degrees while the altitude is in meters. The amplitude is in V/m 
and the time is in hours: minutes on July 1 1, 2002. 




Figure 4. Electric field mill outputs during the roll/pitch maneuvers in Figure 3. There were six 
mills on the Altus aircraft designated ‘Mid-Up’ (MU), ‘Mid-Starboard’ (MS), ‘Mid-Down’ 
(MD), ‘Mid-Port’ (MP), ‘Front-Up’ (FU), and ‘Aft-Up’ (AU). The amplitude is in Volts per 
meter (V/m) and the time is in hours: minutes on July 11, 2002. 
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Figure 5. Ideal fields (g x (t), g y (t), g 7 (*), g Q (t)) in the aircraft frame of reference during the 
roll/pitch maneuvers from Figure 3. The g Q (t) component is estimated using the mean values of 
the 6 mill outputs. The amplitude is in V/m and the time is in hours: minutes on July 1 1, 2002. 
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21:45 21:46 21:47 21:48 21:49 21:50 


Figure 6. Ideal fields (top three plots) and mill outputs (bottom 6 plots) during a period when the 
Citation aircraft was not making calibration maneuvers. Note the features in the waveforms that 
are not due to roll and pitch changes. If we attempted to calibrate the aircraft mills using this 
section of data, the resultant matrix would have almost nothing to do with the true matrix 
M . Amplitudes are in V/m and the time scale is hours: minutes. 
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Figure 7. Ideal (g x (t),g Y (t),g z (t),g Q (t)), gray waveform) vs. realistic ( f x {t),f Y {t),f z (t),f Q {t)) 

, black waveform) fields during a flight of the North Dakota Citation. The amplitude is in V/m 
and the time is in hours: minutes on June 25, 2001. The black plots are the realistic fields based 
on the actual FM outputs while the gray plots are the ideal fields based only on the Gish field 
profile and the aircraft pitch and roll maneuvers (offset by 10 V/m to facilitate comparison). 
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